Tunneling and orthogonality catastrophe in the topological mechanism of 

superconductivity 
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We compute the angular dependence of the order parameter and tunneling amplitude in a model 
exhibiting topological superconductivity and sketch its derivation as a model of a doped Mott in- 
sulator. We show that ground states differing by an odd number of particles are orthogonal and 
the order parameter is in the d-representation, although the gap in the electronic spectrum has no 
nodes. We also develop an operator algebra, that allows one to compute off-diagonal correlation 
functions. 



1. In the BCS theory of superconductivity numer- 
ous physical quantities like London's penetration depth, 
the gap in the electronic spectrum, the tunneling am- 
plitude, etc. are expressed through a single object - an 
off-diagonal two-particle matrix element between ground 
states with N and N + 2 particles 



A(r - r') 



(N\c a (r)c a ,(r')\N + 2) 



(1) 



This fact is not special to superconductivity, but rather 
a manifestation of the mean-field character of BCS. In 
electronic liquids where the interaction is strong in the 
high temperature regime above T c , one expects also to 
see a difference between these different implementations 
of superconductivity. 

Below we consider the extreme case of a strongly in- 
teracting system where a ground state (and the entire 
spectrum) depends crucially on the number of electrons 
in the system. In recent years it has been argued that 
under rather general assumptions, an electronic system 
where a topological soliton is adiabatically attached to a 
particle develops superconductivity. Realizations of this 
phenomenon in low dimensions are Frohlich ideal con- 
ductivity in ID [|] and anyon superconductivity in 2D 
||. We refer to this phenomenon as topological super- 
conductivity . The crucial features of the mechanism 
are (i) an electron acquires a geometrical phase in field of 
a soliton, and a related phenomenon: (ii) orthogonality 
catastrophe. Both of them make the physics of the super- 
conducting state drastically different from BCS physics 
and in particular give rise to an angle dependence of the 
matrix element ([!]) and the tunneling amplitude. 

In this Letter we consider an ideal 2D model of a doped 
Mott insulator on a square lattice at doping close to the 
half filling. In this model the Fermi surface consists of 
four pockets around kf = (±7r/2, ±tt/2). We calculate 
the Josephson tunneling amplitude and show that the 
phase difference between the tunneling amplitude on the 
different faces (1,0) and (0,1) of the crystal is ir. At the 
same time the gap has no nodes. This result is depicted 
in Fig.l and is expressed by eq. (p^). Although it is in 
agreement with the corner-SQUID-junction experiment 
pj the order parameter is different from the conventional 
d-wave form. In a forthcoming paper Q we will show that 
at incommensurate doping the phase difference between 



order parameters in the points (1,2) and (3,4) is the twice 
the angle 8 between vectors ri^ and ^4. 

2. Let us start with a general comment about Joseph- 
son tunneling in the presence of orthogonality catas- 
trophe. Let one side of the junction be a BCS su- 
perconductor with a small gap and phase </>q. Then 
the standard Ambegaokar-Baratoff formula gives for the 
Josephson current I ~ Ime"* ^ F(p„,w)du/w 
in terms of the F-function of a topological superconduc- 
tor at momentum p n normal to the surface, F{p, ui) — 
2iriJ2(N\c p1 \N + l;p)(N+l;p\c- pl \N + 2)6{w - e p ). In 
BCS theory the main contribution to the integral comes 
from intermediate states with an energy of the order of 
a gap, i.e. a pair is destroyed while tunneling. In a topo- 
logical superconductor the ground states \N) and low 
energy states \N + 1) have different topological charges 
and are orthogonal to each other. Their overlap vanishes 
in the macroscopical system. The F-function remains 
nonzero due to a small contribution of infinitely many 
states with energy much larger than the gap. As a re- 
sult, F decays slower than uj^ 1 (in fact F(p, u>) ~ w~ 1//2 ), 
so that J dpF{p 1 lo) is not very small at large u>. Thus, we 
conclude that the tunneling amplitude is given by the the 
equal time two-particle matrix element (|l]) - in contrast 
to BCS, electron-pairs remain intact while tunneling (for 
a similar phenomenon see ^|). 

3. A number of successive steps j^,[i],[l0| have been 
made over last years towards the derivation of the topo- 
logical mechanism from an electronic model with an in- 
finite on-site repulsion 



7^2. 
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Below we sketch such a derivation and add some new 
features to take into account spin correlations. The t — J 
model does not have distinct scales to isolate the physics 
of topological fluids. To capture the physics of interest 
we employ an adiabatic approximation, i.e. we treat a 



hole's motion in a slowly varying spin background |11|. 

A single hop of a hole destroys the short range anti- 
ferromagnet order. However, two consecutive hops and a 
spin flip bring the antiferromagnet in order. As a result, 
the spin configuration remains approximately unchanged, 
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only in second (even) order in the hopping process. In or- 
der to treat spins adiabatically we must first integrate out 
these virtual processes, so that holes remain on the same 
sublattice (a Schrieffer-Wolf like transformation) jlO],[| . 

In second order of perturbation theory in t the hopping 
Hamiltonian is 

U = ^{4(a)W(a,a')^(a') + 4WW(M')c ff '(6')} 

where the sum runs over neighbors (a, a') of sublattice 
A and (b, b') of B on a square lattice. The variable hop- 
ping amplitudes, say i CT(T '(a, a') depend on the spin con- 
figuration and are superpositions of chiralities over paths 
connecting points a, a': 

W(a,a') = t^W aa/ {a,b,a'), (3) 

b 

W aa ,{a,b,a') = (I+ ff -S )(i+ff-S 6 )(i + <r-S /) 

where the overall scale for hopping is t ~ t 2 /J. Now the 
hopping Hamiltonian is ready for the adiabatic approxi- 
mation. 

First we must find a static spin configuration, i.e. the 
amplitudes t(a, a') to minimize the energy. Their modu- 
lus is determined by the competition between electronic 
and magnetic energies and gives the overall scale of the 
model. As far as the phase of the hopping amplitudes is 
concerned we assume that it is determined by the elec- 
tronic energy alone [ fi"l|| . The flux hypothesis |j) sug- 
gests that, in the leading order in doping, the energy 
achieves its minimum if the chiralities along contours 
a — > a + ej — > a + 2e^ are equal in all directions, while 
chiralities along two different paths connecting sites on 
the diagonal of a crystal cell have a different sign (the 
principal of maximal interference) and cancel the ampli- 
tude for diagonal hopping: 

(W(r,r + e x ,r + e x + e y ) + W(r,r + e y ,r + e x + e y )) = 0, 
(t(r, r + e x + e y )) = 0, (t(r, r + 2e 2 )} = l/m. 

The Fermi-surface of the mean field state consists of 
four pockets around Dirac points k/- = fc±.± = (±5, ±5), 
so we decompose electron operator onto four smooth 
movers c CT (r) = 2~2 kf c o--±±( r ) e4k/r - I n what follows we 
refer to the smooth functions c a .±±(r) as to the contin- 
uum part and to the factors e zk/r as to the lattice part of 
the fermion operator c CT (r). 

In this basis, the mean field Hamiltonian can be writ- 
ten in a continuum limit as the square of Dirac oper- 
ator H = D 2 = (a x id x + a y id y ) 2 where the 4x4 
Dirac matrices {a x ,a y } — act in the space labeled 
by (±±). The choice of these matrices (gauge free- 
dom) corresponds to a relabeling of the Dirac points 
and is limited by the symmetry group of the Fermi sur- 
face (there are only four different gauges). We choose 
them to be a x = t% Cg) r^,a y = t\ <E> T3 , where the first 
Pauli matrix r acts on the first (x) ± label and the 
second acts on the y ± label. They correspond to the 



Landau gauge on the lattice with hopping amplitudes 
a x = a(r,r + e^) = l,a y = a(r,r + e v ) = (~l) x . 

Now we are ready to take into account smooth fluctu- 
ations of the phase of the hopping amplitudes (fluctua- 
tions of moduli are not that important). The most effec- 
tive way to do this is to introduce a non-gauge-invariant 
operator ip a to describe charge motion and U(2) gauge 
field A + A ■ cr to describe fluctuations of spin chirality 
Wrr^i ~ expz(F + T ■ a 1 ), where F and T are fluxes 
of U(l) and SU(2) gauge fields. We find the hopping 
Hamiltonian to be Pauli operator: 

H = ^-^{(iV-A-A-ct) 2 +T(F + F-(t)}^ (4) 

The second term with T — —ia x x a y describes the di- 
agonal hopping due to the fluctuations of chirality. 

5. The perturbative vacuum (where the gauge field A 
is small), is unstable when we start to dope the system. 
The energy achieves its minimum if the Abelian part F 
of the flux (the topological charge of magnetic solitons) 
is equal to the number of dopants. The reason for this 
is that the non negative hopping Hamiltonian (|J) in the 
presence of a static flux acquires states with zero energy 
Jl2| . The number of zero modes is twice the number of 
flux quanta and the density of particles occupying zero 
modes is adiabatically coupled to the flux: p(r) = 2^. 
The wave functions (non gauge invariant) of zero modes 
are $(r) = e l ^~ Tx g where A z = d z {—i<t> + x) an d g (the 
lattice part of the zero mode) obeys Tg = —g. If the flux 
is directed up there are two solutions 

9A(B) = e lk++ ' r ± e lk — r + l ( e ik -+' r ± e lk +- r ) (5) 

which are chosen, such that gA = on sublattice B and 
gB = on sublattice A. 

While doping, electrons want to create and occupy zero 
mode states to minimize their energy. This effect com- 
petes with the magnetic energy of the flux. Below we 
assume that the gain in electronic energy wins this com- 
petition . As a result two electrons with opposite spins 
may occupy the same zero mode state. Once zero mode 
states are occupied, the interaction between them lifts 
the degeneracy, so that the zero mode states form a nar- 
row band. In a singlet state the band is always completely 
filled and is detached from the rest of the spectrum, so 
that the chemical potential lies in a gap Ao- This re- 
sults in superconductivity — a density modulation can 
propagate together with a flux configuration without dis- 
sipation, while the electronic spectrum has a gap at the 
Fermi surface. 

A short range antiferromagnetic interaction would sug- 
gest that electrons with spin up (down) spend most of 
their time on sublattice A(B), so that for low energy 
states, we may identify the spin and the sublattice and 
to project onto zero mode states. Thus we obtain the 
U(l) x U(l) anyon model 
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H = —^W — A — A 3 a 3 ) 2 }^ + 2A — + -{T 3 ) 2 
2m In l\ 

[A x (r),A y (r)]=i^5(r-r'). (6) 

The last term with a phenomenological constant A is 
added to implement the magnetic part of the model. 
It induces an attraction between particles with opposite 
spins and determines the gap Ao and T c . 

Here we do not calculate the structure of this narrow 
band but rather concentrate on the distribution of elec- 
tronic spin within the band. It is governed by the SU (2) 
part of the flux. In the "unitary" gauge, where T is di- 
rected along the third axis T — T 3 , the spin density is 

p 3 (r) = l/2ct a 3 c = T 3 /2tt. (7) 

An important supplement to this theory is the identifica- 
tion of the original electronic operator c a on the lattice 
and the non gauge invariant tj) a in the continuum. The 
gauge invariant electron operator creates a hole plus the 
flux of the Abelian gauge field attached to it, i.e. it con- 
sists of a product of ip„ and the vertex operator V, which 
unwinds the Abelian gauge field: H(A)V = (V*)" 1 ?^). 
It is V(z) — exp{— J r A z (z')dz'}, where z = x + iy is a 
holomorphic coordinate relative to the crystal axes. The 
vertex operator creates flux [V(r), F(r')] = irV{r)8(r — 
r'). Thus c a (r) — V^(r) , cr (r)C/(r)(; (J (r) where the fac- 
tor U(r, C) = Ylc a i i s the lattice part of the Dirac tail 
along some contour. It makes the wave function of the 
zero mode gauge invariant. The Hamiltonian (^|) together 
with the correspondence between continuum and lattice 
fields is the field theory for a doped Mott insulator on a 
bipartite lattice. 

6. We now proceed with matrix elements. Let us add 
two particles in a singlet state into a state with N par- 
ticles. The ground state with N + 2 particles consists 
of two extra electrons in the zero mode state and also a 
proper redistribution of the flux T 3 : 

|JV + 2)=£ CTCT , J drdr'fi(r,r')cUr)cUr')\N)*(r,r') (8) 

Operator fi(r, r') creates the flux of the A 3 field according 
to (|): [T 3 (u),fj,(r,r')] = TrfJ,(r,r')(6(r -u)- 6{r' - u)), 
and , 3 / (r, r') is the wave function of a singlet in the zero 
mode state. A solution for the flux operator is /z(r, r') = 
^3( r )^3~ ( r ') wnere Vs(z) is the vertex operator 

V 3 (z) ~ e~ /" A3 *( z "> dz ' = J ln t—^-psiz^dz'dz' ^ 

and L is the size of the system. The meaning of this 
result is simple. The form of operator /x suggests that 
in the spin singlet state electrons with spin a see +27T 
flux attached to the other electrons with the same spin 
while —2tt attached to particles with the opposite spin 
|l3f . This can be illustrated by the operator algebra: 



V 3 (z)c hl (z>) ~ (^)±V* Cu (z')V 3 (z) (10) 
By means of ([h]) we obtain 

A(*, z') = *(z, z')(N\»(z, z'))\N). (11) 

The last factor here is merely a constant while two others 
contribute to the angle dependence. 

Let us first note that the size of the system dropped 
out and is replaced by a short distance cutoff a. This 
shows that the overlap of the ground states with N and 
N + 2 particles does not vanish in a macroscopic system 
|p^| . In contrast, an attempt to insert a single electron 
into the system leaves a zero mode unfilled and creates 
a non-singlet excitation. This leads to an orthogonality 
catastrophe: Due to the operator algebra (|l0| ) the matrix 
element (N\ C<7 \N+l) - (N\c a {0) J dzV a {z)cl{z)\N) - f 
vanishes as L — > 0. 

7. The two particle wave function ^(r,! - ') in (||) con- 
sists of a lattice part and the smooth BCS wave function 
of the pair *(r, r') = */ attice (r, r')A BCS (r - r') 

A SC s(r) ~ / dke 1 ^ A ° 

J v/(e k - /i) 2 + A£ 

where k is relative to minima k±± of the mean field spec- 
trum efe, and A is a gap which separates the narrow 
band from the spectrum. The lattice part of the wave- 
function is ^i a in ce 

(r,r') = U{v,C r )U{v',C' r )g A {v)g B {v') 
depends on two strings (contours) ended in points r 
and r'. Fluctuations of the string are physical exci- 
tations of the pair (not an artifact of the approach). 
In the commensurate case, strings fall in four groups 
within which ^ \ a m ce is the same. These groups cor- 
respond to the states with pairing from different Fermi 
points ky and k^, i.e. to a pair with a total momen- 
tum P = k/ + k^= (0,0) = (±7r,±7r),(±7r,0),(0,±7r). 
A particular string of the wave function with momentum 
P = can be chosen as two contours following each other 
from some reference point up to the point r = (x,y) and 
then a single string along the y-axis to {x,y') and finally 
to the point r' = (x',y') along the x-axis. In the cho- 
sen gauge this factor is U(r, C r )U{r', C' r ) = {-\fiv-v'). 
Then the order parameter (^) is translational invariant. 
Combining this and (|^) we obtain: 

where R = r r' = (X, Y). The numerator of this 
expression is a discrete analog of the continuous holo- 
morphic function in the denominator. Under 7r/2 rota- 
tion it produces the factor e~ l7r / 2 . Another e~ l7r / 2 factor 
is produced by the continuum part. Both phases add 
to e _47r = —1, so that the tunneling amplitude belongs 
to an irreducible d-representation of the crystal group 
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A(-Y,X) = -A(X,Y), It is instructive to look at the 
tunneling amplitude in momentum space. It is: 

A(k) = ^e- iarg ( k ^- targ ( k ~ k ^/(|k-k / |), (13) 

where arg(k/) = j, ^ and f(p) is a smooth func- 

tion. The tunneling amplitude consists of two vortices - 
one in the center of the Brillouin zone (lattice part) while 
another is at a Fermi point (Fig.[|). 
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FIG. 1. The order parameter in momentum space: A(k) 
consists of four similarly oriented unit vortices around each 
Fermi point (±5, ±5). The phase of each vortex is relative 
to the direction of the k/ as in eq. ([TJ|). It changes the sign 
under 90° rotation A(—k y ,k x ) — — A(k x ,k y ). 

8. Eq.(p^[) suggests an interesting generalization to 
the incommensurate case (optimal doping) where the 
Fermi surface is a simply connected curve rather than 
four Fermi pockets. We rewrite ([TJ|) approximately as 

e — iarg(r) r 

A(r) = j- : J dfce lk - r e- larg ( k )A BCS (k) 

Here k is relative to the center of the Brillouin zone and 
the factor arg(k±±) is replaced by arg(k). It is possible, 
because AbcsO*) se * the integral to the Fermi surface. 
At small doping it is shown in Fig.|l|. Away from half- 
filling 6^ — [i ~ i'/(|k| — kf). Thus we obtain a "tomo- 
graphic" representation of the order parameter [ jl5[ : 

A(k) = e- i2aTS ^ J dk f A BCS (k -k f - d)D k ^ f (q)dd 

where the propagator D p (q) = | q | 1 exp (ipq) is a holo- 
morphic function of q. Here pq is an angle between p 
and q and the integral dkf goes over the Fermi surface. 

This is the main result of this paper. It clarifies the 
physics of a topological superconductor. In contrast to 
BCS, an electron with momentum k close to kj emits 
soft modes of density modulation with the propagator 
Dk-k f (q). As a result: (i) the ground states differing by 
an odd number of particles are orthogonal; (ii) the BCS 
wave function is dressed by soft density modes. This is 
analogous to ID physics and bremsstrahlung of QED. The 
new features are: (i) emission of the soft mode is forward; 



(ii) the phase of the matrix element of the soft mode is 
the angle relative to the Fermi momentum. 

9. Although the order parameter ( |i~3| ) forms a 
d-representation, its form and the physics behind it 
are drastically different from the conventional d-wave 
cos k x — cos k y . Nevertheless, it seems premature to spec- 
ulate on an observable difference, until intcrlaycr tunnel- 
ing is taken into account. This is crucial, since the sign of 
parity breaking alternates between odd and even layers 
p6|,p[, and a realistic junction averages over many lay- 
ers. In a hypothetical monolayer tri-crystal experiment 
one would expect the trapped flux to be an integer (in 
contrast to half- integer for conventional d-wave UW). 
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